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Abstract 
Oh! 
D . 

• We compute bounce solutions describing false vacuum decay in a $ model 

in two dimensions in the Hartree approximation, thus going beyond the usual 
^ i one-loop corrections to the decay rate. We use zero energy mode functions 

^ I of the fluctuation operator for the numerical computation of the functional 

determinant and the Green's function. We thus avoid the necessity of dis- 
cretizing the spectrum, as it is necessary when one uses numerical techniques 
based on eigenfunctions. Regularization is performed in analogy of standard 
perturbation theory; the renormalization of the Hartree approximation is 
based on the two-particle point-irreducible (2PPI) scheme. The iteration to- 
wards the self-consistent solution is found to converge for some range of the 
parameters. Within this range we find the corrections to the leading one-loop 
approximation to be relatively small, not exceeding one order of magnitude 
in the total transition rate. 
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1 Introduction 



One of the possible mechanisms at work in the cosmology of the early uni- 
verse is the transitions between two different ground states, a metastable one 
in which the universe may be trapped, and a stable one, the true vacuum 
Pl El El in E] As in an infinite volume system a vacuum tunnelling or thermal 
over-the-barrier transition cannot occur globally, such processes occur by the 
spontaneous formation of regions of "true vacuum" inside the metastable 
state. The transitions may occur either by quantum tunnelling or by ther- 
mal over-the-barrier transitions. If they occur by vacuum tunnelling in a 
space of infinite spatial extension, the classical solution which describes the 
local tunnelling process is called "the bounce" . For a D-dimensional space- 
time the solution is a solution in D-dimensional Euclidean space, which is 
SO{D) symmetric. Here we will consider (1 -|- l)-dimensional space time, 
and, therefore, a bounce in 2 Euclidean dimensions. 

In leading order the transition amplitude is given by jHl 13 IHl El 

r^-'°°P = ^V-'/'expi-S,,) = |^exp(-Sei - S,^,) . (1.1) 

where and ^i-i are the classical and one-loop effective actions, respec- 
tively. The prefactor arises from the functional integration over the transla- 
tion mode riQ = N''V(j), where (p is the classical solution. By a virial theorem 
the normalization of this mode is related to the classical action via 

= Sel . (1.2) 

The computation of the classical action and of the one-loop quantum 
corrections has been performed by several groups for various models |im ITT| 
IT^ IT^ IT^ [T3] The quantum action is of course divergent. Its renormaliza- 
tion should and can be done in exact analogy to the standard renormalized 
perturbation theory. This has been been emphasized in particular in Refs. 
[ini Cni UZl UH Uni where in this, or related, contexts the Born series expan- 
sion has been used to separated the divergent parts in a Lorentz-covariant 
way. 

One may ask whether this leading-order formula, which is based on the 
semiclassical expansion, will receive strong corrections in higher orders. This 
will of course depend on the precise system and on the parameters. For 
electroweak bubble nucleation, a thermal over-the-barrier transition, these 
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corrections were found to be strong ^21; for the parameter sets which at the 
time were beheved to be reahstic, and next-to-leading order corrections have 
been considered by Siirig [201 ■ What he computed were bubble solutions 
which minimized not the classical but the entire one-loop effective action. 
This is different from other calculations [21] in which the bubble was cal- 
culated using the effective potential and quantum-corrected kinetic terms. 
Siirig's calculations were performed for a coupled channel problem involving 
Higgs and gauge fields, which somewhat conceals the simplicity and technical 
elegance of the approach, demonstrating on the other hand that it can be 
used for such involved and realistic problems. More recently the problem 
of computing self-consistently corrections beyond one-loop has been taken 
up by Bergner and Bettencourt |22] who use the Cornwall- Jackiw-Tomboulis 
(CJT)|23j or two-particle-irreducible (2PI) approach. They include the first 
nonleading loop correction, which results in the Hartree approximation. The 
same approximation also occurs as the lowest nontrivial loop correction in 
the so-called two-particle-point-irreducible (2PPI) approach of Coppens and 
Verschelde [211 123 ■ In comparison the Siirig's approach, the Hartree approx- 
imation takes into account the back reaction of the quanta not only to the 
classical field, but also to the quantum fluctuations themselves. Bergner and 
Bettencourt have considered such self-consistent corrections for two systems: 
the kink one space dimension and the bounce solutions in two Euclidean di- 
mensions. Here we will reconsider their approach, but use a different method 
for computing the quantum corrections, analogous to Siirig's, but generalized 
to the Hartree approximation. Both the Green's functions and the fluctua- 
tion determinant will be computed using mode functions of the fluctuation 
operator at zero frequeny, thus avoiding a cumbersome summation over eigen- 
modes. Unlike Bergner and Bettencourt we find, within our parameter range, 
that the fluctuation operator has an unstable mode and a mode at almost 
zero frequency which can be identified as translation mode, both of which 
then can be dealt with in analogy to the one- loop case [TUl 121] • For the 
parameter values of Ref. [22] our iterative computation of the self-consistent 
configuration does not converge, so that a direct comparison is not possible. 

While the computation of the corrections to vacuum tunnelling in the 
Hartree approximation is new in Ref. [22j the Hartree approximation has 
been used to a great extent in finding selfconsistent fermionic lumps describ- 
ing the nucleon in the chiral quark model P71 ISB]. There techniques similar 
to the ones discussed here can and have been applied [29 , 30 

The plan of the paper is as follows: in section [21 we recall the basic 
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formulae of the classical bounce solution. In section |21 we formulate the 
Hartree approximation, on the basis of the 2PPI formalism. In section 0] we 
discuss the translation mode problem. In section |S1 we explain our way of 
computing the Green's function, in section IHl we describe the computation 
of the fluctuation determinant. The numerical aspects of the unstable and 
translation modes are presented in section [7| Renormalization is discussed 
in section IHl Our numerical results are the subject of sectional Conclusions 
and an outlook are presented in section [TUl 

2 The bounce 

Let us consider a scalar field theory in 2D, with the Lagrange density 



and displays two minima, one at $ = 0, and the other one at $ = $_ > 0. 
The value of the potential at this second minimum is lower and represents 
the "true vacuum" , while $ = represents the "false vacuum" . 

The bounce is an SO (2) symmetric classical solution of the Euclidean 
field equations [t —>■ —it). We denote the Euclidean variables hy xi = x and 
X2 = T. The radius is r = \/ xf + x^, and the classical field is denoted as 
(f){r). The classical Euclidean action is given by 




where the potential U ($) is given by 




(2.2) 




(2.3) 



and the bounce which minimizes this action satisfies 



A20 + f/'(0)=O, 



(2.4) 



or 



(f) ld(j) 



+ t/'(0) = 



(2.5) 



r dr 
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and boundary conditions 



^|,=o = 0, 0™ = $+. (2.6) 
dr 

The one-loop correction to the classical action is given by 

where m is the mass in the false vacuum, 

= f/"(0) . (2.8) 

and where the prime denotes that the translation zero mode is removed 
and that one replaces the imaginary frequency of the unstable mode by its 
absolute value. 

The transition rate from the false to the true vacuum is given, in the 
one-loop approximation, by 

r^-'°°P = V-^/^ exp(-5ei) = (1^) exp(-5,i - 5i_i) . (2.9) 

The prefactor arises originally as the normalization of the zero mode. The 
presence of a zero mode in this approximation is demonstrated by taking the 
gradient of the classical equation of motion: 

Vi [-A20 + U'm = [-A2 + f/"(0)] = . (2.10) 

Its normalization, defined by rj^i = AfoVicj), and the condition that rjoi is 
normalized to unity, is given by 

Afo' = J d'x{S/4f , (2.11) 

where there is no summation over i. The right hand side is, for a spherically 
symmetric solution 0(r), equal to the kinetic term. Furthermore, one finds 
that in two dimensions, using a scaling argument, that J d'^xU{(j)) = when 
is the classical solution. We therefore obtain Af^'^ = Sd, i.e., Eq. (|1.2p . 
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3 The bounce in the Hartree approximation 

The Hartree approximation can be derived from the 2PPI formahsm as fol- 
lows: The effective action in this formalism is given by 

S,4M', <P] = SM + r'^'^^'lM', - Y / d'xA\x) , (3.1) 

up to renormalization counter terms discussed in section 8. Here A is a local 
insertion into the propagator which has the form 

g'\x) = -A2 + M\x) , (3.2) 

with the definition 

=m^ + ^X(f)^ -6r](f)+^XA . (3.3) 
iA itself is defined by the gap equation 

Finally r^^^^ is the sum of all two particle point irreducible graphs, in which 
all internal propagators have the effective masses Ai"^. A graph is two particle 
point reducible (2PPR) if it falls apart if two lines meeting at a point are cut. 
To lowest order in a loop expansion r^^^^ is given by a simple loop, i.e. 

r-" = ll,.det4i±4. (3.5) 

2 -A2 + m2 ^ ^ 

We will come to the zero mode question later. In this approximation, A is 
given by 

= =< ''' -A, + A^2 l^^>= ^(^^^) ■ (3-6) 

Here the Green's function Q is defined by 

i-A2 + M')gix,y)=5'ix-y). (3.7) 

In taking variational derivatives of the effective action we have to consider 
A as a function of Ai^ and 0, i.e., in the last term of Eq. ()3.H1 we have to 
replace 

A = -4>^ + ^ [M^ - + 6r]^) , (3.8) 

OA 



byEq. 

Taking the variational derivative of the effective action with respect to 
the field then leads to 




(3.9) 



Using rotational symmetry we obtain explicitly 




x|=r 



The back reaction of the quantum modes onto themselves is contained in 



In the formula for the transition rate the functional determinant was modified 
by removing the zero mode, whose functional integration is not Gaussian 
and which leads to the prefactor Sci/2n. Once we modify the fluctuation 
operator by introducing the quantum back reaction, there will be no zero 
mode anymore. However, if Q{x,y) were the exact Green's function 



the presence of the zero mode would still follow from the classical equation 
of motion STes/Scplx) = and translation invariance. As an artefact of the 
approximation the zero mode is shifted to some finite value ^ 0. This 
is a well known problem, which also arises in other applications. So in the 
Hartree approximation to the 0(4) sigma model, the pions do not have mass 
zero, in spite of their role as Goldstone bosons jHH IS21 tHB] As long as the 
corrections are small, the lowest eigenvalue of the fluctuation operator in 
the partial wave I = 1 will still be close to zero, as found by Surig in his 
calculations on bubble nucleation. Of course, if we do not remove this mode, 
the Green's function will receive a huge contribution from it. Furthermore, 
removing a factor of dimension (energy)^ from the functional determinant 



4 The translation mode 




(4.1) 
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is necessary in order to provide the transition probability with the correct 
dimension (energy)"^. 

We take here a pragmatic point of view in removing the "almost zero" 
mode from the determinant. The formalism then requires to remove the 
translation mode contribution from the Green's function as well, as A is now 
the functional derivative of the modified functional determinant. 

Of course identifying and removing the "would be" translation mode in 
this way constitutes an approximation which comes in addition to the Hartree 
approximation itself. We will not be able to trust this approach if either of 
these approximations leads to large modifications of the bounce and of its 
determinant. In the case of the zero mode we will expect to satisfy 
\uj^\ << AE where AE is the typical level spacing. As most of the modes 
are continuum modes, the "typical" level spacing would be the energy differ- 
ence between the zero mode and the unstable mode. We will continue this 
discussion in section 9 in conjunction to our numerical results. The technical 
problem of removing the pole will be dealt with in the section 8. 

A more fundamental approach has been taken long ago, when the mass 
of the sine-Gordon soliton received great attention. Various schemes |34| 
ESI IMl EH EH EH] were formulated to deal with the problem of translation 
invariance and the projection to zero momentum states. These have been 
used in order to compute the two-loop correction to the soliton mass [101 • 
They become very involved in higher loop orders, and we are not aware of 
a formulation for the case of a resummation as the one considered here. We 
therefore have not followed such an approach. 

As mentioned above, the prefactor 5'ci/27r originates from integrating out 
the zero mode, and by using virial theorem relating the normalization of the 
zero mode to the classical action. In the Hartree approximation this virial 
theorem no longer holds. Nevertheless from what we have discussed above, 
the exact zero mode should be V0ci, where 0ci now is the self-consistent 
profile. We therefore will use its normalization as the prefactor, so that the 
transition rate in the Hartree approximation becomes 

pHartree ^ If^j^^^SeS ^ (4.2) 

where the effective action is defined in Eq. (j3.1|) . 
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5 Computation of the Green's Function 



In order to include the back-reaction of the quantum fluctuations to the 
bounce in the Hartree approximation we need the Green's function Q{x,x') 
of the (new) fluctuation operator. In fact the Green's function is usually 
discussed in a more general form, as a function of energy. Here such a concept 
corresponds to introducing an additional third dimension. We will choose it 
spacelike, thus introducing an Euclidean time. Using translation invariance 
in the time direction we introduce the Fourier transform Q{x, a;', i/^), where u 
is the Euclidean frequency. The generalization, equivalent of introducing an 
additional dimension, is necessary for discussing the translation mode and 
reappears in the formulation of the determinant theorem in the next section. 
The Green's function satisfies 

[-A2 + m^ + V{r) + u^]g{x, x') = 5\x - x') , (5.1) 

with 

V{r) = -6r](^{r) + ^\{(P\r) + g{x,x)) . (5.2) 

The Green's function can be expressed by the eigenfunctions of the fluctua- 
tion operator. We denote them by ipa{x), they satisfy 

[-A2 + m^ + V{r)]M^) = ulilJaix) . (5.3) 

In terms of these functions the Green's can be written as 

^, 2^ ^a(^)^a(^') 

We may, furthermore, decompose the Hilbert space into angular momentum 
subspaces, introducing eigenfunctions exp{ilip)Rni{r), where (p is the polar 
angle and where the radial wave functions Rniir) are eigenfunctions of the 
partial wave fluctuation operator: 



Then the Green's function takes the form 



Rni{r) ^ uliRniir) . (5.5) 



^(f,f' ) = Y^Y^,ai^-^') Rni{r)Rni{r') ^^^^ 
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These expressions are formal, we have discrete and continuum states, so the 
sum includes summation over discrete states and integration over contin- 
uum states. While these expressions arc very suitable for discussions on the 
formal level, they are not very suitable for numerical computation. In partic- 
ular, if one uses these expressions for the numerical computation, it becomes 
necessary to discretize the continuum states by introducing a finite spatial 
boundary. 

There is a well known alternative way of expressing Green's functions. 
Consider first the free Green's function obtained for V{r) = 0. It can be 
written as 

and this may be expanded as 

^ oo 



(5.7) 



(5.8) 



l=—oo 



where r< = min|a;|, r> = max|x|, \x'\ and = + z/^. Note that 
ultimately uo"^ will be zero and k = m, so we have to deal with the modified 
Bessel functions /; and Ki. They satisfy 



1 d? 



/2 



Blinr) = , 



(5.9) 



where Bi stands for Ii or Ki. Ii{K,r) is regular at r = and Ki^nr) is 
exponentially decreasing for r — > oo. Their Wronskian is given by 

Ki{Kr)dIi{Kr) / dr — Ii{tir)dKi{Kr) / dr = 1/r . 

We now expand the exact Green's function in an analogous way by the ansatz 

^ oo 



(5.10) 



l=—oo 



The functions f^{r, v"^) satisfy the mode equations 



+ ^ + K + V{r) 



dr"^ r dr'^ 



/z^(r,^^) = 0. 



(5.11) 
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and, furthermore, the following boundary conditions: 



fi (r, z^^) oc r' r ^ ,^ 

fj^(r,iy'^) (X exp{—Kr)/y/l^ r ^ oo ' 

So f[~ is regular at r = and /^^ is regular, i.e., bounded, as r — >• oo. For 
V{r) = these boundary conditions are those satisfied by Ii{nr) and Ki{Kr), 
respectively. Furthermore, as the behaviour at r = is determined by the 
centrifugal barrier, and the behaviour for r cxd by the mass term, these 
boundary conditions are independent of the potential. If we write 

/r(r) = hiKr)[l + h;ir,u')], (5.13) 
/+(r) = KiiKr)[l + hiir,u')], (5.14) 

then the functions /if (r, u'^) become constant as r — and as r — >■ oo, and 
for finite r they interpolate smoothly between these asymptotic constants. 
If we impose the boundary conditions /i^(r, i/^) — >■ for r — »• oo then the 
Wronskian of fj^ and f[~ becomes identical to the one between Ki{Kr) and 
hi^r), i.e., equal to 1/r. Applying the fluctuation operator to our ansatz, 
Eq. ()5.1()|1 . we then find 

11 °° 1 

[-As + + V(r)+] G(x, x', z/2) = —-6(r-r') V e^'^^"^') = -6(r-r')6(^-^') . 
L J 2tt r ^-^ r 

l = — OD 

(5.15) 

This completes the construction of the Green's function. 

Numerically we proceed as follows: the functions hf" satisfy 

which can be solved numerically. The second differential equation is solved 
starting at large r = f with h'l{f,i>'^) = hi (f, i/^) = 0, and running back- 
ward, f has to be chosen far outside the range of the potential. In this region 
/if(r, z/^) are essentially constant. For the first differential equation we first 
obtain a solution hiir^p'^) starting at r = 0, with /^^(O, i/^) = /ij(0,z/^) = 
0. This function does not satisfy the boundary condition required for the 
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Green's function, it will be used for the computation of the functional deter- 
minant. The function /i^~(r, i/^) is obtained from hilr,!/^) via 

K(ry ) = '2 f r . (5.18) 

l + hi{r,iy^) 

which obviously solves the differential equation with the appropriate bound- 
ary conditions. 

Finally the Green's function is given by 

^ oo 



=— oo 



(5.19) 



6 Computation of the Fluctuation Determi- 
nant 



The fluctuation determinant which appears in the rate formula 



V = det'- 



-A2 -|- m? ' 



(6.20) 



can be written formally as an infinite product of eigenvalues of the fluctuation 
operator. The prime denotes taking the absolute value and removing the 
translation mode. As in the previous section we introduce the generalization 

P(.=,.dat4i±4±^. (6.21) 

Note that we omit the prime, here. Using the decomposition of the Hilbert 
space into angular momentum subspaces we can write 



)=n 



Ln 



UJ 



ln{0) 



n 

1=0 



det M,(z/2) 
detMj°V) 



di 



(6.22) 



with the radial fluctuation operators 



dv^ r ar 



(6.23) 
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as before, di denotes the degeneracy. If we restrict / to positive values then 
di = 2 ioi I > and di = for I = 0. 

According to a theorem on functional determinants of ordinary differential 
operators ^ |32| we can express the ratios of the partial wave determinants 
via 



det Miii^^) 



det Mf°)(z/2) 



lim 



where il)i{v'^,r) and %l)f'\v'^ ,r) are solutions to equations 

mi{v^)Uv\r) = , mf\v^)^ljf\v\r) = . 
with identical regular boundary conditions at r = 0. Of course 

Furthermore we have 



Iiinr) 



(6.24) 



(6.25) 



(6.26) 



(6.27) 



where hi{v'^,r) is precisely the function we have introduced in the previous 
section, except for the fact that we have added to the fluctuation operator. 
We finally have 

detMi(z/2) r . 9 ^ . X 



det m[°V) 



1 + hAv^ , oo) 



and 



InPfz/'l 



1=0 



1 + hAv^ , oo) 



(6.29) 



The fluctuation determinant in the transition rate formula and in the 2PPI 
formalism refers to the fluctuation operators at z/^ = 0, and so in the nu- 
merical computation we just need the functions /^^(O, oo), as for the Green's 
function. The only exception is the translation mode we will discuss in the 
next section. 
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7 Unstable and translation modes 



In the one-loop formula for the transition rate the determinant of the fluc- 
tuation operator appears as det'(— A2 + A1^), and the prime denotes two 
modifications with respect to the naive determinant: 

(i) the unstable mode has an imaginary frequency, corresponding to a 
negative eigenvalue ul = — z/^ of the fluctuation operator. It is to replaced 
by its absolute value. This mode appears in the s-wave I = and manifests 
itself by a negative value of l + ho{0, 00). So here we have to take the absolute 
value. 

(ii) the translation mode manifests itself, in the one-loop approximation, 
by the asymptotic limit 1 -|- hi{—uj^,oo) = 0. We denote the frequency of 
the translation mode, which is the lowest radial mode in the m = 1 partial 
wave, by uiq = ut. The fluctuation determinant ()6.21|) has a factor uj^ + 
which has to be removed, according to the definition of det'. Otherwise the 
logarithm of this expression, appearing in the functional determinant, does 
not exist. Furthermore the Green's function is not defined either at = —lv^. 
In the one-loop approximation = 0, in the Hartree approximation it is 
close to zero, otherwise we cannot trust our approximation of identifying this 
mode with a "would-be" zero mode. 

In the one-loop approximation the translation mode is removed numeri- 
cally in the following wav|lUj: we compute /;.i(oo,±e^) for some sufficiently 
small e and replace 



i.e., we take the numerical derivative at u"^ = 0. 

As discussed above, beyond the one-loop approximation there is no exact 
translation mode, but a pole in the Green's function appears very close to 
1/^ = 0, and its wave function is close to V0. This contribution would 
make the corrections extremely large. As long as the corrections beyond 
one-loop are small, the "almost zero" mode still corresponds to a collective 
motion of the system, with almost vanishing restoring force. As discussed in 
section 3 we will continue to remove it also in the Hartree approximation. 
We determine the position of the eigenvalue by requiring 1 + hi{—uf, 00) to 
vanish, and compute the numerical derivative not at = but at u"^ = —ujf, 
i.e., we remove a factor uj"} + . 



l + /ii(0, 00) 




(7.1) 



2e2 



14 



The Green's function in the / = 1 channel has, at r = r', the form 



We can use the fact that the pole term is antisymmetric with respect to 
+ iJl by computing the Green's function at = —oj"^ ± and by taking 
the average of these two values. Then the pole term has disappeared and the 
averaged Green's function takes the form 

1 , ,2 _ , ,2 

\r> ^^ ^ , .2 , ,2\ I _ , .2 ,2\1 d2 /^\ ^In 



I !2 , ,2\2 ^4 

(7.3) 

As long as uo^ and are much smaller than the tu^^ this is a good approxi- 
mation to the desired reduced Green's function 

[^i(^^0)U = J]^. (7.4) 



In 



In the explicit numerical computation of the Green's function we use of 
course the expression ()5.10|) . As evident from Eqs. ()5.18|) the pole arises 
from the re-normalization of the mode function /ii(z/^,r), i.e., from dividing 
by 1 -|- h_{i>'^, oo). In averaging over the Green's functions at z/^ = —uj^ ± 
we add two very large terms which almost cancel. This can be done in 
a somewhat smoother way: if is sufficiently small we can assume that 
1 -\- hilu"^ , oo) passes through zero linearly and we may replace 

1 + hi{-ujt ± , oo) ^ ±^ hi{-ujt + , oo) - hi{-uj^ - , oo) . (7.5) 
The average over the Green's functions can then be cast into the form 
. , n^l ~ /i+(-c^? + e\ r)U-^I + ^^ r) - ft{-uj - e\ r)h{-uj^ - 

[yi{',',^)\red — r / 2,2 \ I f 2 2 \ 

ni[—uj^ + e"^, oo) — ni[—ujf — e^, oo) 

(7-6) 

where /i(z/^,r) = /i(/tr)[l + hi{i''^,r)] is the mode function before the 
re-normalization. 
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8 Renormalization 



In the previous sections we have presented the basic formahsm and its nu- 
merical implementation. There is still one point to be discussed: divergences 
and renormalization. It is well known, that renormalization in (1 + 1) dimen- 
sions just requires normal ordering. This means in practice that we have to 
redefine A — Q{x, x) by 

A{x)^g{x,x)-go{x,x), (8.7) 
or, in terms of the functions hf{r) 

oo 

A(f) = diIi{mr)Ki{mr) [hl{r) + hj{r) + hl{r)hj{r)\ . (8.8) 

1=0 

Normal ordering is not unique, it depends on the mass used in the free 
Green's function. We have used here the mass in the false vacuum; this may 
be changed, but it essentially means that we redefine our couplings r) and A. 

There is one more place where we have to subtract a tadpole diagram: in 
the tr log term. If we expand it perturbatively we find 

d^xV{r)gQ{x,x) . (8.9) 

While it is trivial to remove this part in a perturbative calculation it is 
less obvious how to remove it from a nonperturbative one. Indeed we have 
computed In V using the partial waves. However the divergent contribution 
can be traced exactly in this partial wave representation: it is given by the 
contributions of first order in V{r), and these can be computed exactly. 
We have 

+ [2m|^ + -]^}Hr) = V{r)[l + k{r)] • (8.10) 
ar"^ Ii[mr) r dr 

This inhomogcneous differential equation allows for an iterative expansion of 









-InV 






2 


div 





hi (r) and the first order contribution h\ (r) is the solution of 



(i^ I'Amr) 1, d 



0-2 + P^yTH + -]l-yhf\r) = V{r) . (8.11) 
dr^ Ii[mr) r dr ^ 
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This equation can easily be solved numerically. Then the renormalized con- 
tribution of a partial wave / to the fluctuation determinant is given by 



Ji = di 



+ hi{oo)) - hf\oo) . (8.12) 



This renormalization can be generalized to higher dimensions, using the 
renormalization procedure of Verschelde |S1 which applies to the Hartree 
approximation. Then higher subtractions are required, they can be likewise 
determined exactly within the numerical procedure, see, e.g., ^21 El]- While 
Verschelde' s procedure of removing the divergences is based on an elaborate 
analysis of Feynman graph's, it can be, from the practical point of view [HH] , 
obtained by adding a counterterm 

S^,t. = j d^x [AM\x) + A] . (8.13) 

A is the "cosmo logical constant" counterterm, in (3 + 1) dimensions there is 
a further counterterm proportional to Ai^{x). Finiteness of the gap equation 
and of the action require 



2 , 47r/i 



A = --7 + ln^ , (8.14) 

87r [e J 

A = -m^A . (8.15) 



This corresponds precisely to the subtractions we have made, which are iden- 
tical to minimal subtraction. In particular 

A{x) = g{x,x) + A = &mte . (8.16) 

The expressions for the Green's function and for the one-loop action imply 
summation over angular momentum I. After renormalization these converge 
in the expected form, the single terms being proportional to while the 
unsubtracted terms behave like This has been verifled numerically and 
presents a cross check on the procedure and the accuracy. We have appended 
the sum from L = /max to oo by adding this sum in its asymptotic form. We 
have used a flt A/l^ + B /l^ + C/l^ through the flve values of J/ at Zmax — 
5, . . . , /max in order to parameterize the terms in this asymptotic sum. In 
order to illustrate the convergence in / we plot, in Fig. Q the values Hl and 
of where the latter is the exact flnite sum up to L complemented by 
the asymptotic sum, for A = 0.9, rj = 0.5. The convergence is seen to be 
excellent. In our actual calculations we used /max = 25. The same procedure 
was applied to the summation over angular momenta occuring in G{x, x). 
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Figure 1 




Figure 1: Convergence of the angular-momentum summation, a = 0.9,(3 = 1; 
dots: S^; dashed hne with crosses: 

9 Numerical results 

The computational framework introduced in the previous sections has been 
carried through for a representative part of the parameter space. The pa- 
rameter space is conveniently analyzed ITO] in terms of two dimensionless 
variables a and (3 which are introduced as follows: We define the scaled 
variables: X = mx and $ = 27] /m?(j). Then the classical action takes the 
form 



4 

m 



= /5^ei(a), (9.1) 

with (3 = /4:7f and a = \(3/m?. The variable a can take values < a < 1 
as a condition for the existence of the classical bounce solution, a = 1 
corresponds to degenerate minima, the limit a — > 1 is called the thin-wall 
limit. This defines the parameter range for a. (3 appears in front of the 
rescaled classical action while there is no such factor in front of the quantum 
action which in the one- loop approximation only depends on a. Therefore 
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large (3 imply a large classical action and therefore relatively small quantum 
corrections, small values of (5 make the classical action small and the quantum 
corrections relatively important. This consideration applies to the one- loop 
approximation; for a self-consistent scheme the separation of classical and 
quantum parts is not unique, but this consideration may still serve as an 
estimate. Indeed we find that for (3 < 0.8 the Hartree iteration, even with 
an underrelaxation parameter, does not converge. 

The computations were started with the one-loop approximation and then 
iterated until the largest difference in the profiles 0(r) between two subse- 
quent iterations max^ A0(r) was smaller than 10~^. The convergence of the 
profiles 0(r) and of the potential V{r) is displayed in Figs. |21^a) and (b), 
respectively, for the case a = 0.9,(3 = 1. 




Figure 2: Convergence of the bounce profiles and of the potential for a = 0.9 
and P = 1. 

The iteration is found to converge for a < 0.9 and f3 > 0.8. The transla- 
tion zero mode is found at typically |z/^| ~ 10~^ in the one- loop approxima- 
tion, thus verifying to good accuracy the quality of the classical profile and 
of the integration of the mode equation. In the Hartree approximation the 
mode which we continue to identify with the zero mode is located at values 
of |z/^| ~ 10~^, the largest value of 0.1 occurs for a = 0.2,(3 = 0.8. For 
a = 0.9 and /? = 1 the required accuracy is obtained after 80 iterations, this 
takes about 40 seconds with a 1.3 GHz processor. For a < 0.9 and (3 > 1 
the convergence is much faster. 

The a-dependence of the unstable mode in the one-loop approximation 
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is displayed in Fig. 01 One sees that it approaches zero for « — > 1. This is 
obviously the cause of the lack of convergence for a > 0.9. During the iter- 
ation towards the self-consistent Hartree solution, the mode moves between 
positive and negative values of u"^ while giving large contributions to tr^. If 
we look, see Fig. |31 at the effective actions found in the Hartree approxima- 
tion then there is no evidence of any singular behaviour near a = 0.9. So we 
think that the lack of convergence is just a technical problem which could 
be overcome, but this certainly would need some new analytic idea, not just 
plain numeric efforts. 




Figure 3: The unstable mode for [3=1. Straight line: one- loop approxima- 
tion; squares: Hartree approximation. 

The dependence on the parameter a at fixed (3 is displayed in Fig. 0] 
The classical actions and the quantum actions in the one-loop and in the 
Hartree approximations are plotted including a factor (1 — a). Indeed all of 
them are singular in the thin- wall limit a 1. For the classical action in the 
one-loop approximation the behavior can be determined analytically using 
the standard technique (see e.g. as 

lim5r°°^(l-a)=/5-. (9.2) 

The one loop results extend to a = 0.97, for the Hartree results we were not 
able to obtain convergence beyond a > 0.9 as already mentioned. 
In Figs. Efa-c) we plot the quantum actions 
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Figure 4: Classical and total actions as functions of a, multiplied by (l — a). 
Dashed line: 5'^^"'°°''; empty square: thin wall limit of 5'^^"'°°''; straight line: 
^to7^°°^; full squares: S^^'''''^''; full circles: 5^^^'''". 



as obtained in the one-loop approximation and 

obtained in the Hartree approximation as functions of p. It is found that the 
one-loop action, which is independent of (3, and 5'^^''*'''^'' approach each other 
for large values of (3. When (3 decreases below (3 = 1 the difference increases 
rapidly, pointing, for a = 0.2 and 0.5, to a kind of singularity at values of 
f3 between 0.8 and 0.9. Below this range the iteration ceases to converge. 
During the iteration the value of 0(0) decreases, and the classical profile 
displays a minimum at some finite value of r. The lack of convergence for 
P < 0.8 is not well-understood. For some parameter sets, e.g., a = 0.5, Fig. 
Efb), there seems to be a kind of singularity of the quantum action, which, 
however, is compensated by an opposite behavior in the classical action, so 
that is not apparent in their sum, as can be seen in FigE^a). Some analytic 
idea, like reshuffling contributions between the quantum and classical parts, 
could possibly cure the problem. Indeed in a self-consistent scheme there is 
no real "classical" part, as the profile depends on the quantum corrections. 

In Figs. infa,b) we display the difference between the one-loop and the 
Hartree approximations for various quantities, for a = 0.5 and a = 0.9. 
ASq = 5'^'^''*^^'^ — 5'q~^°°P, the difference of the quantum actions, is found to be 
negative for a = 0.5 (and also for a = 0.2), but positive for a = 0.9. We also 
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Figure 5: Quantum parts of the effective action as functions of (3. Dasfied 
lines: one-loop approximation; squares: Hartree approximation; (a): a = 0.2; 
(b): a = 0.5; (c): a = 0.9. 
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Figure 6: Differences between one-loop and Hartree approximation of various 
parts of tlie effective action as functions of 13. Notations as specified in tlie 
text; (a): a = 0.5; (b): a = 0.9. 
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display the difference between the classical actions ASd = S^'^^^'^^'^ — Sl^^°"^ 
and the total difference AS'tot = AS'd + ASq. We see that even this total 
difference has different signs for a = 0.5 and a = 0.9. We should keep in 
mind again that the separation between classical and quantum action is not 
unique and that the classical profile, and therefore the prefactor, depends on 
the quantum fluctuations. The total width, which includes in the prefactor 
the normalization of the zero mode is for all a found to be suppressed in the 
Hartree approximation with respect to the one-loop approximation. This is 
also displayed in Figs. Efa) and (b), where the solid line represents 

pHartree 

^logr = -log^^, (9.5) 

with the one-loop transition rate given by Eq. (jl.lll . and rate in the Hartree 
approximation given by Eq. ()4.2j) . This difference is positive for both a = 0.5 
and a = 0.9. These rates are displayed separately in Fig. [7| 

We have mentioned above that for the transition rates we use the nor- 
malization of the translation mode as prefactor, as the virial theorem used 
in the one-loop approximation is no longer valid in the Hartree approxima- 
tion. In the one-loop approximation the normalization of the zero mode and 
the classical action agree better than four significant digits, in the Hartree 
approximation the classical and total actions approach the normalization of 
the zero mode for large f3, but /3 ~ 1 both of them differ from it by factors 
up to 3. 

As we have mentioned above the procedure becomes unstable for (3^1; 
the iteration ceases to converge. Bergner and Bettencourt 22j, on the other 
hand present results for m = 1.3469, and, in our conventions A = 2 and 
rj = 0.9685. They use the same renormalization convention as we do. Their 
parameters correspond to a = 0.967 and /3 = 0.877. This is close to the 
thin wall limit, and at a value of (3 for which our procedure to becomes 
unreliable. The features of the solution found in Ref. [S] are quite different 
from those found here: there the eigenvalue spectrum displays neither a 
zero nor an unstable mode. It may be that the tunnelling in this region, 
f3 < 0.9, a > 0.9 takes place in a qualitatively different way. If this is so then 
our procedure certainly breaks down on more than numerical reasons as it is 
based on an conventional treatment of zero and unstable modes. In order to 



^The potential of is given by V{<j)) = \[(jP - v'^ f/A - e{(f> + v)/{2v). In order to 
compare it with ours it has to be shifted in (j) such that the left minimum is at (/> = 0. 
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le-15l [ • ^ • J 

Figure 7: The transition rates in the one-loop and Hartree approximations 
as functions of f3. One-loop approximation: straight line: a = 0.2; dashed 
line: a = 0.5; dash-dotted line: a = 0.9; Hartree results: circles: a = 0.2; 
squares: a = 0.5; diamonds: a = 0.9. 

obtain a comparison of the result of Ref. [22^ with ours, we have determined 
the various actions for (3 = 0.877 and m = 1.3469 fixed, as a function of 
a. These results are plotted in Fig. |H1 The classical action of Ref.. [22] 
agrees reasonably well with our results. We also show the thin wall limit, 
Eq. (|9.2p . The total one-loop action can be computed up to a = 0.98. The 
total one-loop action obtained in Ref. |22] is somewhat lower than our values. 
Surprisingly our results for the Hartree approximation seem to extrapolate 
quite naturally towards the result of Ref. j22] • 

10 Conclusions and outlook 

We have presented here a calculation of the transition rate for the false vac- 
uum decay in (1-1-1) dimensions in the Hartree approximation. We have 
used analytic and numerical methods that have been applied previously to 
other systems, for the computations of functional determinants [101 HH I2E] 
and zero point energies fHl 121] ■ Both techniques are based on mode func- 
tions of the fluctuation operator at z/^ = 0, as they appear in the standard 
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Figure 8: Comparison with Ref. |22]: Dashed hne with stars: our results for 
Sl^^°°^; full line with stars: our results for Sl^^°°^; empty circles: our results 
for S^^, full circles: the corresponding results of Ref. [221; empty square: 
thin wall limit of the classical action; all actions are multiplied with (1 — a). 

representation of Green's functions. Here these techniques are used for the 
first time in conjunction; as one sees this leads to a rather simple and effec- 
tive computation scheme with various cross checks. We should like to add, 
that these techniques are well-suited for coupled-channel systems as well [12] . 
Renormalization can be dealt with in exact correspondence to standard per- 
turbation theory; this feature extends to applications in higher space-time 
dimensions and, allowing for dimensional regularization, to computations in 
gauge theories. 

We have presented results for a substantial part of parameter space, using 
the parameters j3 = m^/Arf' and a = \l3/m?. Bounces exist only for a < 1, 
a = 1 being the thin wall limit. The second parameter 13 weighs the ratio 
between classical and quantum contributions to the effective action. We find 
that for /? > 0.8, and for a < 0.9 our iterative determination of the field 
configuration in the Hartree approximation converges. The unstable mode 
and an almost-zero mode in the p-wave, identified as the translation mode, 
persist during the iteration and are handled as in the one-loop approximation. 
We find that the transition rates are generally suppressed in the Hartree 
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approximation, as compared to the one-loop rates. However the corrections 
are relatively small, less than an order of magnitude. This has been found 
similarly for the one-loop quantum corrections for instanton transitions in 
the two-dimensional Abelian Higgs model |^ and may be an general feature 
of two-dimensional models. 

The only similar computation which is available at present is the one of 
Ref. [22|- Unfortunately our procedure does not converge for the parameters 
chosen by these authors, which in our convention are a = 0.967 and /3 = 
0.877. Their self-consistent configuration is qualitatively different from ours 
in that the fluctuation operator has no unstable and translation modes. If 
this is the case in this part of the parameter space then it means that the 
tunneling proceeds in a qualitatively different way, there. Nevertheless we 
have found that our Hartree results seem to extrapolate naturally towards 
the ones of Ref. [221 ■ 

The methods used here for the computation of self-consistent solutions 
for bounces naturally extends to bounces in higher space dimensions and to 
other classical solutions in quantum field theory. In particular, the techniques 
for covariant regularization and renormalization are available PSl QH] in the 
framework we have established here. 
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